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ABSTRACT 

We find a strong correlation between the peak energy at zero fluence (E pea k,o) an d the isotropic 
energy (E 7iso ) of the 22 pulses of 9 Gamma Ray Bursts (GRB) detected by the Fermi satellite. The 
correlation holds for the individual pulses of each GRB, which shows the reality of the correlation. 
The derived correlation (Spearman correlation coefficient, r, is 0.96) is much stronger compared to the 
correlations using E pca k (in place of E poa k,o) determined from the time-integrated spectrum (r = 0.8), 
or the time-resolved spectrum not accounting for broad pulse structures (r = 0.37), or the pulse-wise 
spectrum (r = 0.89). Though the improvement in the E poa k - E 7 .j so relation (the Amati relation) for 
a pulse-wise analysis is known earlier, this is the first time a parameter derived from a joint spectral 
and timing fit to the data is shown to improve the correlation. We suggest that E pca k,o, rather than 
E pca k, is intrinsic to a GRB pulse and a natural choice as the parameter in the pulse- wise correlation 
studies. 

Subject headings: Gamma-ray burst: general — Methods: data analysis — Methods: observational 



1. INTRODUCTION 

Gamma Ray Bursts (GRBs) are the most luminous 
events in the universe, predominantly observed in the 
hard X-ray to gamma ray energies. An individual GRB 
appears as a flash of gamma ray event lasting for seconds 
and then continuously shifts towards the lower energy 
bands - all the way to radio wavelength (see for example, 
Meszaros 2006). In the initial phase of high energy emis- 
sion, known as the prompt emission, it undergoes sig- 
nificant spectral evolution. Over the years, researchers 
have developed empirical models to describe the time- 
integrated spectrum, the light curve and the spectral evo- 
lution. It has been shown that various model parameters 
correlate with the energy related physical parameters of 
the full GRB. For example, the parameter peak energy 
of the v¥ v spectra (E pca k) correlates with the isotropic 
energy (E 7jiso ) (Amati et al. 2002), isotropic peak lu- 
minosity (Lj SO ) (Schaefer 2003; Yonetoku et al. 2004), 
and collimation-corrected energy (E 7 ) (Ghirlanda et al. 
2004). Similar correlations hold in the time domain e.g., 
spectral delay {n a g) - L; so (Norris, Marani k, Bonnell 
2000), variability (V) - L; so (Fenimore & Ramirez-Ruiz 
2000) and rise time (r rise ) - L iso (Schaefer 2007). The 
objective of such correlation studies is two-fold. Firstly, 
to use GRB as a standard candle of cosmology, alongside 
the commonly used standard candle - supernova. Sec- 
ondly, to understand the GRB physics itself. GRB can 
be observed at a much higher redshift (as high as z > 6) 
than supernova (z < 1.7). With the advent of the Fermi 
satellite we now have the record holder - GRB 090423 
(z = 8.2; see Ghirlanda et al. 2010). In spite of this 
tremendous advantage over the supernova, GRB has one 
serious drawback. Unlike the supernova, whose mecha- 
nism is very well understood and supported by numer- 
ous observations, the physics of GRB is not understood 
properly. Hence, it is interesting to investigate empirical 
models which can predict the correct energetics of GRB. 



The predictions of theoretical model should conform to 
that of a data driven empirical model and this whole pro- 
cess gives a strong constraint on the possible physics of 
GRB. 

Pulse analysis of GRB prompt emission has received 
considerable attention in recent times. It has the poten- 
tial to unravel physical processes responsible for the ob- 
served correlations and help in standardizing the energy 
budget so that GRBs can be used as precise cosmolog- 
ical distance indicators (Hakkila 2009). In the prompt 
emission phase, a GRB varies vastly in time and en- 
ergy domains, so much so that no two bursts have the 
same temporal and spectral characteristics (Norris et al. 
1996). This situation simplifies if the whole GRB event is 
considered as an ensemble of temporally separated pulses 
of self similar shapes (Nemiroff 2000), generated nearly 
simultaneously in a wide range of energy bands (Norris 
et al. 2005). Hence, each of the pulses can be mod- 
elled separately and then the whole GRB event can be 
reproduced by shifting and adding these pulses. One of 
the strongest constraint the pulse analysis gives is that 
the pulses of a GRB, despite having separate set of pa- 
rameters and different energetics, have the same redshift. 
Hence, each pulse can be used as a distance indicator and 
must conform with each other. It has been shown that 
all the correlations studied for the full GRB event hold 
good, sometimes even better, if the inherent pulse prop- 
erty of GRB is taken into account. Krimm et al. (2009), 
for example, showed that the E pea k - E 7! ; so correlation 
derived from pulse-wise analysis is consistent with the 
time-integrated analysis. Hakkila et al. (2008) studied 
Tiag - Lj so correlation and conclusively showed that the 
spectral lags are pulse properties rather than burst prop- 
erties. 

These kind of correlations are, however, empirical in 
nature, and may not have any physical significance (Band 
& Preece 2005; Butler et al. 2007; 2009). Moreover, they 
do not use the full information available in the sense that 



2 



the spectral and timing parameters are considered inde- 
pendent of each other. In such correlation studies, the 
information of time evolution of the spectral parameters 
and the energy dependence of the timing parameters are 
lost. Hence, the derived spectral and timing parameters 
are average quantities. Liang & Kargatis (1996) showed 
that one of the spectral parameters, namely E pea k, fol- 
lows a time evolution law. Kocevski & Liang (2003) used 
this interdependence to study the evolution of E pca k, and 
found that E pea k exponentially decreases with the fluence 
of an individual pulse in a GRB. Recently, Basak & Rao 
(2012) have developed a new method for the complete 
empirical description of the individual pulses of GRB 
prompt emission, simultaneously in the time and energy 
domains, based on well established empirical formulae. 
The method conclusively shows that the two new pa- 
rameters of the Liang & Kargatis (1996) model, namely, 
the peak energy at zero fluence (E pe ak,o) and the char- 
acteristic fluence (</>o), can be used to predict the timing 
parameters, namely, spectral delay and width. In other 
words, these two parameters are more intrinsic to the 
GRB pulses than the average spectral and timing pa- 
rameters, and the other pulse characteristics can be de- 
rived from them. In this paper we demonstrate that the 
pulse-wise analysis gives a better correlation than both 
the time-integrated analysis and time-resolved analysis, 
which do not take broad pulses into account. We also 
show that Ep 0a k,o, instead of E poa k, gives even better 
correlation, and it is the correct choice of parameter in 
correlation studies. The structure of the paper is as fol- 
lows. In §2 we briefly describe the method and in §3 we 
give the data analysis and results. Major conclusions are 
discussed in §4. 

2. SIMULTANEOUS TIMING AND SPECTRAL 
DESCRIPTION OF A GRB PULSE 

Pulse analysis of GRB essentially involves extracting 
the constituent substructures based on various empirical 
time description of pulses (Norris et al. 1996; Nemiroff 
et al. 2000; Hakkila et al. 2009; 2011), e.g., the fast rise 
exponential decay model (FRED; Kocevski 2003), the 
Norris model (Norris et al. 2005). Though this method 
is useful in extracting isolated and slightly overlapping 
pulses, it cannot be used to extract "heavily-overlapping" 
and low signal-to-noise pulses at present (Hakkila et al. 
2011). The energy spectrum of a GRB is popularly 
described by an empirical model given by Band et al. 
(1993). The same model can also be used for the indi- 
vidual pulses. Choosing Norris model for the time do- 
main and the Band model for the energy domain, and 
employing the time evolution of the peak energy (E pea k) 
of Band model, as proposed by Liang & Kargatis (1996), 
Basak & Rao (2012) developed a technique to determine 
the model parameters (E poa k.o and <p ). 

Basak & Rao (2012) analyzed the Swift/BAT and 
Fermi/GBM data of the brightest GRB in the Fermi era 
- GRB 090618. They measured the global parameters of 
the individual pulses of this GRB, namely, the low energy 
photon index (a) , high energy photon index (/?) using the 
Band model (Band et al. 1993) for the time-integrated 
spectral data, and the characteristic time scales (ti and 
T2) by fitting the energy- integrated light curve with the 
Norris model (Norris et al. 2005). These global parame- 
ters were then used to generate the XSPEC table model 



with the parameters Epea^o and <fto as variables. From 
spectral fitting in XSPEC, E pe ak,o and <j>o were deter- 
mined. 

The data analysis for the global parameters (see Basak 
& Rao 2012) showed improvement due to the inclusion of 
the Swift/BAT along with the Fermi/GBM. The simulta- 
neous spectral fitting, however, showed some systematic 
errors in the overlapping energy regions. Hence, they 
used only Swift/BAT data for the simultaneous timing 
and spectral description of the GRB pulses. The typical 
value of the peak energy (E pea k) of GRBs (and the indi- 
vidual pulses) is ~ 300 keV, and we expect still higher 
values of the parameter, E pea k.o- The Swift/BAT energy 
range (15 - 150 keV; Barthelmy et al. 2005) is inade- 
quate for accurate measurement of these parameters in 
many occasions. Hence, the derived values by Basak and 
Rao (2012) showed large uncertainties (typically 50 - 100 
keV). GBM onboard the Fermi satellite provides an un- 
precedented energy range (~8 keV to ~40 MeV) with 
adequate sensitivity (Meegan 2009; Meegan et al. 2009) 
for constraining the values of E pca k and E poa k,o- Hence, 
in the present analysis, we have used Fermi/GBM data. 
Also, we have carefully chosen the time interval of each 
pulse avoiding contamination from the other pulses. For 
example, in GRB 090618 (see Table 1), we have chosen 
time intervals 61-76, 76-95 and 106-126 seconds for the 
pulses 2, 3 and 4 respectively (compared with the time 
divisions of Basak & Rao (2012): 50-77, 77-100 and 100- 
130 seconds, respectively). We have included the precur- 
sor burst of this GRB (pulse 1) in the present analysis. 

3. DATA ANALYSIS AND RESULTS 

In our analysis, we have chosen the set of GRBs which 
were considered for time-resolved spectral analysis by 
Ghirlanda et al. (2010). This contains 12 long GRBs 
with known redshift, detected by Fermi/GBM till the end 
of July 2009. Among these, three GRBs (GRB 080905, 
GRB 080928, and GRB 081007) are very weak bursts 
and could be fit only with a single power law with un- 
constrained peak energy. Hence, we have used the data 
for the remaining 9 GRBs and their individual pulses (a 
total of 22 pulses). Table 1 contains the full list of our 
sample with the name, measured redshift, and pulses of 
GRBs in the first, second, and the third columns, respec- 
tively. 

In a given GRB, we select those pulses which have 
broader width compared to the rapid spikes in its light 
curve. While selecting such pulses we carefully avoid 
those pulses or portions of the pulses which have too 
much overlap with others. Fishman et al. (1995), for ex- 
ample, have discussed various categories of light curves: 
1. single pulse, 2. smooth, cither single or multiple, 
well-defined peaks, 3. distinct, well-separated episodes 
of emission, 4. very erratic, chaotic and spiky burst. 
More than one such category of pulses can show up in 
a single GRB. In such cases, the fourth category gives 
rise to many overlapping pulses in some parts, while the 
other parts of the same GRB might be dominated by 
the second or third category of pulses. The temporal re- 
gions populated by the fourth category do not allow a 
unique measurement of model parameters of a pulse, if 
each pulse has to have an independent set of parameters. 
Hence, we avoid such regions and take only clear portions 
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TABLE 1 

The observer frame peak energy (E pea k), zero fluence peak energy (-E pea fc,o) AND the isotropic energy (E Jz i 

INDIVIDUAL PULSES OF GRBs (GRB DATA SET TAKEN FROM GHIRLANDA ET AL. 2010) 
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of a burst. 

The broad pulses selected in our analysis are listed in 
Table 1 (also see Ghirlanda et al. 2010: Figure 2), along 
with the appropriate start time (t\) and stop time (£2) 
shown in columns 4 and 5. Ghirlanda et al. (2010) have 
done a time-resolved analysis, whereas we have done a 
pulse-wise analysis. In the former case, the GRB light 
curve is arbitrarily divided into a large number of bins 
and the evolution of spectra is examined. Krimm et al. 
(2009) have shown that the correlation between the peak 
energy and the isotropic energy of broad GRB pulses im- 
proves from that of the whole GRBs. In our analysis, we 
follow the same approach and select the time-cut accord- 
ing to the broad pulse structure. Some portions of the 
light curves are neglected, because, either these are dom- 
inated by the fourth category of pulses or have low count 
rates to accurately determine model parameters. GRB 
080810 and GRB 080916C have low count rates after 30 
s and 55 s respectively. GRB 090323 contains multiple 
overlapping spikes between 30 to 59 s and 75 to 135 s 
region which do not contain any broad pulse structure 
(also count rate is low in 75-135 s region). GRB 090328 
has two overlapping spikes between 20 to 26 s (taken as 
20-24 s and 24-26 s in Ghirlanda et al. 2010) and hence 
neglected. 

We essentially use the method described in Basak & 
Rao (2012). The global parameters (a, (3, n and T2) are 
determined to generate three dimensional pulse model 
for a set of E poa k,o and t/V The time- integrated spectra 
for these set of values give a two-parameter XSPEC table 
model. We perform x 2 minimization of the spectral fit of 
the data with this model to determine the best fit values 
of the model parameters and their nominal 90% confi- 
dence level errors (A% 2 = 2.7). Also, the normalization, 
which is a free parameter in the model, is determined by 



XSPEC. 

The results of our spectral analysis are shown in Ta- 
ble 1. The observer frame peak energy (E pea k) calcu- 
lated from the Band model, the zero fluence peak energy 
(E P eak,o) an d the measured nominal 90% confidence level 
errors for each are shown in columns 6 and 8. In order 
to compare the improvement in correlation of these pa- 
rameters with E 7i i so (column 10 of Table 1), we perform 
a linear fit of the form log (y) = K + 5 log (x) using the 
technique of joint likelihood for the coefficient K, 6 and 
the intrinsic scatter ((Tint) (D'Agostini 2005; Wang et al. 
2011). Here, y = E poa k (full and pulse-wise study) and 
E P oak,o, are in the units of 100 keV, and x = E 7! ; so is in 
the units of 10 52 ergs. 

In Figure 1, we have plotted the source peak energy 
(E P eak) as a function of E 7i i so for the 22 individual pulses 
of the 9 GRBs (stars). For comparison, in the same 
figure we have plotted the time-integrated values for 
the full GRB (filled boxes), based on the data given 
in Ghirlanda et al. (2010). For illustration purposes, 
a scatter plot of the values for time-resolved spectral 
analysis is also given in the figure (small circles). Note 
that Li SO is the more appropriate parameter for a time- 
resolved spectral analysis, but to compare with our re- 
sults we have converted the flux given by Ghirlanda et 
al. (2010) to E 7j i so using the multiplication factor of the 
time bin. The straight line shown in the figure is the 
log (Ep Ca k) = K + S log (E 7i i so ) fit to the pulse-wise data 
using the joint likelihood method. The time-integrated 
data is also fit by the same technique. The Spearman 
correlation coefficient (r) , probability that the correlation 
occurred by chance (P), the parameters for the linear fit 
and the intrinsic scatter in the data are given in Table 
2. For comparison, the values reported by Ghirlanda et 
al. (2010) for the time-integrated values for 10 GRBs 
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Fig. 1. — Time-integrated, time-resolved (not accounting for 
broad pulse structure) and pulse-wise E VRa h as a function of 
isotropic energy (B 7j j so ) for the 9 Fermi-GRBs with measured 
rcdshift. The large filled boxes denote the time-integrated points 
while small circles denote the time-resolved points (GRBs arc cho- 
sen from Ghirlanda et al. 2010). The stars denote the pulse- wise 
Epeak as determined by the Band model (present work). 

are also given in the table. A comparison of correlation 
and P for the time- integrated (0.80 and 9.60 xlO -3 , re- 
spectively) and pulse-wise (0.89 and 2.95xl0~ 8 , respec- 
tively) analysis shows that there is an improvement in 
the correlation and the reality of the correlation (i.e., 
lower P) in the latter case, favouring pulse-wise anal- 
ysis. Intrinsic scatter of the data ((Tint) per point is 
also reduced (0.225/9 to 0.244/22). This improvement 
in the correlation of E pea k - E 7! ; so for pulse-wise analysis 
is known in earlier works. Krimm et al. (2009), for exam- 
ple, used a sample of Swift-Suzaku GRBs which gave the 
Spearman correlation as 0.74 with a chance probability 
of 7.58xl0~ 5 for a sample of 22 GRBs. This correlation 
improves, for a set of 59 pulses of these GRBs, to 0.80 
with a chance probability of 5.32 x 10 - 14 . 

The time-resolved study not accounting for broad 
pulses shows a very poor correlation (r=0.37; see Ta- 
ble 2). This is expected because, in such time divi- 
sions, broad pulse structure is not considered. Even if 
the spikes are considered as pulses, they are overlapping 
in any such small time division and hence not usable to 
determine E pea k, which can be uniquely associated with 
a pulse. Time divisions taken for a broad pulse, on the 
other hand, has the facility to study the evolution of 
Epeak and uniquely determine the E pO ak,0 m that pulse. 
In our analysis, we ignore these rapidly varying spikes 
and concentrate on the broad pulse structures. 

In Figure 2, we have shown E pea k,o as a function of 
E 7! ; so along with a straight line fit. A comparison of 
Figure 1 with Figure 2 shows immediately that the new 
parameter E pea k,o has a better correlation in the E poa k.o 
- E 7i i so plane compared to the pulse-wise E pea k - E 7 .j so 
analysis, which already shows improvement in terms of 
correlation, P and <T int compared to the time-integrated 
analysis. All these values for the E pea k,o - E 7! ; so fitting 
are shown in the last row of Table 2. The correlation 
coefficient (r = 0.96) is significantly higher than that of 
Epeak - E 7i i so , whether time-integrated (0.80) or time- 
resolved (0.37) or pulse- wise (0.89). The chance prob- 
ability also decreases. As a comparison of correlations 
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Fig. 2. — Pulse-wise -Ej, e afc,0 as a function of isotropic energy 
(Ej,iso) for the 22 pulses of 9 Fermi-GRBs. Data for three indi- 
vidual GRBs (GRB 090328, GRB 090424, and GRB 090618) are 
marked with separate symbols. The solid line is the best linear fit, 
while the dot-dashed lines show the 3cr scatter. 

found in previous works, we quote some of the results 
presented by Krimm et al. (2009): The original Am- 
ati (Amati 2006) catalogue of 39 bursts has r = 0.87 
with chance probability 4.72 xl0~ 13 . A sample of all 
the long bursts (91) shows a poor correlation (r = 0.76) 
with chance probability 4.72 xl0~ 18 . The sample of 59 
GRB pulses of Krimm et al. (2009) has a correlation 
of 0.80 with a chance probability of 5.32 xl0~ 14 . Our 
analysis shows a clear improvement in the correlations, 
whether time-integrated or time-resolved (not account- 
ing for broad pulse) or pulse- wise analysis. 

It is long debated that the correlation such as E pea k - 
E 7i i so might be a result of observational selection effects. 
In particular, the E pea k values are stretched due to the 
multiplication factor (1+z) making the correlation look 
better (Nakar and Piran 2005; Band and Preece 2005; 
Butler et al. 2007, 2009; Shahmoradi and Nemiroff 2009). 
But, there are other researchers who argue in favour of 
the reality of these correlations (Ghirlanda et al. 2008; 
Nava et al. 2008; Krimm et al. 2009; Amati et al. 2009; 
Ghirlanda et al. 2010). If these correlations hold within 
a GRB as for a sample of many GRBs, then we can con- 
clude that the correlations are indeed physical. In Fig- 
ure 2, we have marked the data points for GRB 090424 
(stars) and GRB 090618 (filled triangle) for which we 
analyze the highest number (four) of pulses. The data 
shows that the pulses of the same GRB follow the same 
Epeak, o - E 7i i so correlation. Also, the fact that this cor- 
relation is tighter than E pca k - E 7i i so correlation signifies 
that E pC ak, o is more intrinsic of a GRB pulse than E pca k 
and hence this parameter should be used for such correla- 
tion studies. It is interesting to note that the pre-cursor 
of the GRB 0901618 also follows the same correlation. 
The GRB 090328 is also marked (filled squares). First 
two pulses of this GRB (with higher Epeak^o values) show 
particularly high deviation. An examination of the E pea k 
- E 7i ; so values of these two pulses (Figure 1) also shows a 
much higher deviation. Interestingly, the third pulse of 
this GRB is consistent with the trend. It is known ear- 
lier that the first pulses of GRBs tend to be harder and 
thus deviate from the correlation (Krimm et al. 2009, 
Ghirlanda et al. 2010). But, this deviation is much less 
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TABLE 2 

Statistical analysis for the correlations of E pea k (I to IV) or E pea k.o (V) with E lt i ao and parameters for the linear fit 
for the 9 Fermi GRBs and their individual pulses. The methods are, I: time-integrated study (present work) for 9 GRBs, 

II: TIME-INTEGRATED STUDY FOR 10 GRBs (QUOTED FROM GfflRLANDA ET AL. 2010). Ill: TIME-RESOLVED STUDY NOT ACCOUNTING FOR 
BROAD PULSE STRUCTURES (CALCULATED FROM GHIRLANDA ET AL. 2010). IV: PULSE-WISE STUDY (PRESENT WORK). V: PULSE-WISE 
Epeak.O ~ ^l,iso CORRELATIONTHE SPEARMAN CORRELATION COEFFICIENT (r), CHANCE PROBABILITY (P) AND THE PARAMETERS OF LINEAR 

FIT (K, 8 AND (Ti n t) ARE REPORTED HERE. 

"Method r P K 8 a mt xj e A Aoi ) 



I 0.80 0.0096 0.166±0.080 

II 0.004 0.162±0.085 

III 0.37 0.0095 

IV 0.89 2.95X10- 8 0.289±0.055 

V 0.96 1.60X10- 12 0.640±0.050 



compared to the pulse- wise E pea k in the E pea k - E 7i ; so 
correlation. 

4. DISCUSSION AND CONCLUSIONS 

The strong correlation between the luminosity and the 
parameter obtained from a joint spectral and timing fit 
for the individual pulses of GRBs clearly indicates that 
the basic radiation/ emission process is similar in diverse 
GRBs and the dispersion in the other parameters like 
Lorentz factor, beaming angle etc. are quite minimal. 
An examination of the data in Table 1 and Figure 2 
shows that the dispersion in the correlation for differ- 
ent pulses of a GRB is of the same magnitude as the 
dispersion between GRBs. Hence it is worthwhile to in- 
vestigate the radiation/ emission process in operation by 
making a direct model fit to the data. This would help 
us finding the fundamental parameter/ variable respon- 
sible for the correlation. If such an exercise brings down 
the dispersion in the relation, it may be possible to use 
GRBs as a distance indicator for cosmological purposes. 
It is also interesting to note that the parameters derived 
for the "pre-cursor" in GRB 090618, is consistent with 
the global correlation. 

Asano & Meszaros (2011) have done simulations of the 
spectral and temporal evolution of gamma-ray bursts us- 
ing internal dissipation models and have concluded that 
the models reproduce the Band spectra and also the 
generic time evolution. Dado et al. (2007) have used 
the master formula based on the "cannonball" model 
and have explained the various correlations observed for 
the prompt emission of the GRBs. The present work 
demonstrates that it is possible to fit the data with a 
comprehensive set of formulae describing the temporal 
and spectral evolution of the bursts. Hence, it should be 
possible to directly fit the data with the model predic- 
tions and derive the fundamental quantities responsible 
for the universal correlation. Any model contains basic 
physical assumptions along with other details. A direct 
fit to the data taking individual pulses should segregate 
the basics from the details. For example, in the "cannon- 
ball" model (Dado et al. 2007) the Lorentz factor and 
the viewing angle of a GRB determine most of the prop- 
erties of a GRB pulse. Since the viewing angle would be 
the same for the different pulses of a given GRB, a di- 
rect fit to the spectral and temporal profile can have the 



0.473±0.048 0.225±0.067 0.64 (7) 
0.476±0.079 0.47 (8) 

0.516±0.049 0.244±0.048 0.56 (20) 
0.555±0.050 0.291±0.039 1.04 (20) 

additional constraint of the constancy of this parameter. 

Further improvements for the method described here 
(also see Basak & Rao (2012)) could be made by itera- 
tively including color corrections to the light curve and 
also re-confirming the E pea j< evolution formula. It would 
also be interesting to repeat this exercise for short GRBs 
and X-ray flares, which will give further clues to the emis- 
sion mechanisms responsible for the correlation presented 
here. 

The choice of isotropic energy (E 7i ; so ) over the other 
two physical quantities, namely, the isotropic peak lumi- 
nosity (Lj so ) and collimation-corrected energy (E 7 ) can 
be justified as follows. Peak luminosity is measured based 
on the assumption that the spectral shape at the peak is 
the same as the average shape, which is not very phys- 
ical (Wang et al. 2011). On the other hand, E 7 could 
have been a better choice, as collimation effect is cor- 
rected. Though, in practice the beaming angle (9j) is 
very ill determined (e.g., see the weak constraints in the 
measured Oj by Goldstein et al. 2011). It is evident 
from Table 2 that the intrinsic scatter does not improve 
in the new analysis. This might happen due to some 
inherent assumptions of our model, crucially the hard- 
to-soft evolution and/or the uncertainty in the measured 
rcdshift. A close inspection of Figure 2 reveals that the 
pulses of the same GRB are scattered on the same side of 
the correlation line. This might be due to the fact that 
the correct energy budget is the collimation corrected en- 
ergy and not the isotropic energy. Analysis of 14 pulses 
of 6 GRBs (for which 9j values could be collected) shows 
a correlation of 0.91. We believe that this correlation 
might improve with increasing accuracy of the measured 

To summarize, we have used the method of joint tim- 
ing and spectral description of GRB pulses and found 
that Ep^o is a fundamental parameter in the pulse de- 
scription. 
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